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Abstract 

We present a simple one-dimensional model to describe the evolution of a transverse 
dune field. The model is characterized by the distances between the dunes and 
their heights, which determine the inter-dune sand flux. The model reproduces the 
observation that the dunes in a given field have approximately all the same height. 
We find that this result is independent on the initial configuration of the field, as 
well as on coalescence effects between migrating dunes. The time for the final state 
to be reached is studied as a function of the relevant phenomenological parameters. 
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1 Introduction 



There are over hundred different types of dunes that have been classified by 
geomorphologists according to their shape and size [1] . Dunes can be found in 
deserts, in coastal areas, on the sea-bottom [2] , and even on Mars [3] . Although 
no general theoretical model has been developed up to now to explain this wide 
variety of dune shapes, it is well known that the dune morphology is mainly 
determined by the amount of available sand and by the behavior of the wind 
over the year. The simplest type of dunes occurs when there is not very much 
sand, and when the direction of the blowing wind doesn't change with time. 
These dunes have a crescent-like shape and are called "barchans". Since the 
pioneering work of Bagnold in 1941 [4], several field measurements led to a 
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good understanding of barchan dunes. In particular, it has been observed 
that these dunes move without changing their shape with a velocity which 
is inversely proportional to their height [5]. Fields with hundreds of barchans 
are found for instance in Peru [6,7,8], Namibia [9] and Marroco [10,11]. In a 
recent work [12], a two-dimensional model to describe the collective motion 
of barchan dunes has been presented, which was able to reproduce the same 
kind of spatial correlations between the dunes as in real dune fields. 

As the sand availability becomes larger, parallel dunes perpendicular to the 
wind direction, which are called "transverse dunes" , arise for a wind blowing 
steadily from the same direction. Ensembles of transverse dunes, with crest 
to crest distances ranging from a few meters to over 3 km, are found to cover 
about 40% of all terrestrial sand seas, being common in the Northern Hemi- 
sphere [13,14], and also dominating the sand landscape on Mars [3]. The last 
few years have seen an increasing interest in the investigation of the geomor- 
phological aspects and the evolution of transverse dunes [15,16,17]. In spite 
of this effort, very little is known about these dunes, and several phenomena 
are still not well understood, as for example the effect of coalescence between 
dunes with different velocities [18]. 

This work reports on a simple one-dimensional model for the evolution of 
transverse dunes in a field, which reproduces the observation that these dunes 
have generally the same height and velocity of propagation. This is indepen- 
dent on the initial configuration of the field and on coalescence effects between 
close dunes with different velocities. In our model each dune is considered to 
have zero width and infinite length, but a variable height, which determines 
the associated propagation velocity and that also controls the interaction be- 
tween neighbouring dunes. The paper is organized as follows. In section 2, the 
model is presented; the results for both cases with and without coalescence 
are discussed in section 3; conclusions are made in section 4. 



2 The Model 



We consider a set of N dunes distributed along the x axis, each one with a 
variable height hi, but zero width. The dunes are allowed to move in the x 
direction with a velocity Vi = dxi/dt given by: 



where Xi is the position of the i— th dune, i — 1, 2, . . . , N, and a is a constant 
phenomenological parameter that contains information about the wind speed, 
the grain size, etc. 
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The behavior of the sand flux on a dune field is crucial to the evolution of the 
field dynamics, as interactions between dunes strongly depend on the sand 
transport among them. For a dune that migrates shape invariantly through 
sand avalanches on its slip face, the sand flux over the dune behaves as 
d<fi/dh =constant, where h(x) is the height profile of the dune. For the model 
studied in this paper, we consider the height hi of the dune at the position of 
the crest, and set for the i— th dune: 

</>i = bhi, (2) 

where b is our second phenomenological parameter, again depending on wind 
speed, grain size, etc. Note that equation (2) means that the rate of sand 
transport over the crest is larger for the bigger dunes than for the smaller 
ones. As a consequence, exchange of sand is expected between neighbouring 
dunes with different heights, as a bigger dune tends to loose sand more rapidly 
than a smaller one. 

To model the inter-dune interaction, we first observe that for a two-dimensional 
slice of a dune, the mass of sand in the dune, or equivalently, its area A, in- 
creases with the height of the dune as A ~ h 2 . On the other hand, the amount 
of sand transported depends on the difference between the sand fluxes on the 
crests of the neighbouring dunes, i.e. if dAi is the volume change for the i— th 
dune during a time interval dt, then dAi ~ (dhi) 2 ~ (0j — <f>i-i). Thus, we 
write: 

^ = -« \J\<pi - 0i_i| x sign(0j - fa-i), (3) 

where c is the third parameter of our model that depends on the shape of 
the dune, and the last factor on the right side ensures that more sand is 
transported from the larger dune to the smaller one. As in a real transverse 
dune field, we set the sand influx to point in the positive wind (x) direction 
and assume it to be constant. We call the sand influx (fi , which can be seen 
as an effective dune fixed at the position xq = 0, with a constant height 
ho = (fio/b- Equations (1), (2) and (3) are repeatedly iterated, which gives the 
time evolution of the dune field. 



3 Results and discussion 

We studied the evolution of a set of N dunes with initial positions 2^(0) and 
heights hi(0), for different configurations with N = 10, as shown in Figure 1. 
In the first group of simulations (I), the dunes are equally spaced along the x 
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axis at t — 0, and have the same heights hi(0) = H,i = 1,2, ... ,N (Figure 
1(a)). In the second group (II), we take random values for the heights /^(0) 
chosen from a homogeneous distribution between and 1, but still the dune 
spacing is regular (Figure 1(b)). And in the third group (III), both the values 
for the heights and the positions of the dunes along the x axis are chosen 
randomly (Figure 1(c)). 

The trajectory of the dunes along the x axis and the changes in their heights 
are shown as a function of the time in Figures 2(a) and 2(b) respectively, 
for particular realizations corresponding to each one of the initial states I, II 
and III. The values of the phenomenological parameters are a = 300m 2 /year, 
b = 1500 year and c — 0.05, and the sand influx is chosen to be 0o — 
2.25 x 10 4 m/year « 4.3 cm/min. In the simulations, each iteration corresponds 
to 0.01 year. As can be seen from these figures, after a certain number of 
iterations t ~ T^, all dunes reach the same height /ij(Too) = <p /b = 15m, and 
consequently the same velocity ViiT^) = a/hiT^) = (a x b)/<p = 20m/year, 
independently on the initial configuration. 

In the discussion above, the distances between the dunes were large enough 
that their relative order was the same during the whole evolution of the field. 
However, if a smaller, faster migrating dune is very closely behind a larger, 
slower one, they can reach the same position after some finite time. Dune coa- 
lescence has been studied for instance in Ref. [12], in which a critical distance 
rfmin between two neighbouring dunes was considered, below which they collide 
and become one single dune. In our one- dimensional model, we set d min = 0. 
For simplicity, dune coalescence is implemented here through the emergence 
of a single dune, with height h c = ^Jhf + hf +1 . Consequently, the number of 
dunes N in the field decreases by unity, whenever two dunes coalesce. 

The main plot of Figure 3 shows the variable Xi(t) for a set of iV(0) = 20 
dunes with random values of Xj(0) and /i«(0) for t = 0. The values of O an d of 
the parameters a, b and c are the same as in Figures 2(a) and 2(b), but now 
their distances in the initial state are much smaller. The number of dunes as 
a function of time for the realization in this figure is shown in the upper-left- 
hand-corner inset. Note that this number decreases monotonically until the 
last coalescence between dunes has occured, and the final number of dunes, 
iV(oo) = 3 has been reached. As we can see from the bottom-right-hand-corner 
inset, coalescence of dunes doesn't modify the result that the dunes in the field 
do all have the same height in the final state. 

For a given distribution of dunes in the field at t — 0, the number iV(oo) of 
dunes in the final state tends to be larger for smaller values of the parameter 
a, the other parameters b, c and 0o being fixed. This is because a controls 
the velocity of the dunes, and coalescence is more frequent in a field with 
faster dunes. In the same way, the closer the dunes are in the initial state, 
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Fig. 1. Different realizations of a set of N dunes with positions Xj(0) and heights 
hi(0) at t = 0. In (a) the dunes are equally spaced, and have the same height H 
(group I). Dunes in group II still have regular spacing but random heights (b), while 
in group III, both the positions and heights are random (c). 
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Fig. 2. (a) Trajectory of the N = 10 dunes along the x-axis as function of time, for 
each of the realizations shown in Figure 1. The evolution of the heights of the dunes 
for each case is shown in (b). The values of the parameters used are a = 300 m 2 /year, 
b = 1500 year -1 , c = 0.05, and </>o = 2.25 x 10 4 m/year. See text for details. 

the more they are going to coalesce on average. As an example, we show in 
Figure 4 the behavior of N(oo) as a function of the length L of the system, 
i.e. the initial distance between the first dune and the last one in the field, 
for a set having a fixed number of N(0) = 50 dunes, using the same values of 
the phenomeno logical parameters as in Figure 1. As we can see in this figure, 
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Fig. 3. The main plot shows the time evolution of the position X{(t) for a set of 
N(0) = 20 dunes with coalescence for the values of the parameters used in Fig- 
ure 2. The number of dunes as a function of the time, N(t), is shown in the up- 
per-left-hand-corner inset of this figure. From the bottom-right-hand-corner inset, 
we can see that dune coalescence doesn't change the fact that all dunes have the 
same height for t — > X^. 

the smaller (larger) the value of L, the more (less) the dunes tend to coalesce. 
The statistical error bars are obtained from different realizations of dunes with 
random heights and positions at time t — 0. The full line shown in Figure 4 
corresponds to the fit: 



where N(0) = 50, q ~ 1, and L* ~ 1.34 x 10 3 m. In Figure 5, we show for 
three different values of L the corresponding curves N^t), averaged over five 
realizations of (xj(0), /ij(0)), as a function of time t. The values of L shown 
are L = 100m (circles), L = 500m (triangles) and L = 1000m (squares). 
The dashed line corresponds to a straight line in a linear-log fit to the curve 
corresponding to L = 100m. We can observe that Ni(t) decays approximately 
logaritmically with time, independently on L. After a variable time interval 
T c < T^, coalescence stops and an asymptotic value iV(oo) is reached. 

As all dunes attain the same height ho ~ 0o f° r t = T^, the time to reach the 
final state depends on the strength of the sand flux in the field, <f> . A plot of 
Too as a function of 0o f° r the initial configuration I of Figure la is presented 
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Fig. 4. Average number of dunes N(oo) for t = Too, as function of the size L of the 
system, i.e. the distance between the first dune and the last one in the initial state. 
The statistical error bars are obtained from an average over different realizations 
(xj(0), hi(0)) of N(0) = 50 dunes at t = 0. The full line corresponds to the fit using 
the equation N(oo) = N(0) [1 — <?*exp(— L/T*)], with q ~ 1 and L* ~ 1.34 x 10 3 m. 
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Fig. 5. Three iVx,(t) x t curves corresponding to averages over different initial real- 
izations of the type III (a?j(0), hi(0)) of a set of iV(0) = 50 dunes, which are let to 
evolve with coalescence. The curves correspond to L = 100m (circles) L = 500m 
(triangles) and L = 1000m (squares), where L is the distance between the first dune 
and the last one in the initial state of the field. Time is plotted on a logarithmic 
scale. The dashed straight line is a guide to the eye meant to fit to the curve for 
L = 100m (circles). 
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Fig. 6. Dependence of the time to reach the final state, on the relative sand flux 
\(fro — bH\ for the initial configuration I shown in Figure 1(a). The circles shown in 
the main plot correspond to the case (fro > bH, and line represents equation (5). 
The inset shows as a function of (fro for the same set of dunes. 

in Figure 6. As we can see from this figure, scales after a transient interval 
of variation in (fr — bhi(0): 



i.e. the closer the value of <po to <fi = bhi(0) = bH, the more rapidly the system 
will reach the final state. From equation (3), we see that the phenomenological 
parameter c controls the rate of sand transport between neighbouring dunes. 
We have found that the number of iterations to reach the final state decreases 
as 1/c, which in fact is independent on the initial height H of the dunes. 



4 Conclusions 

We presented a simple one-dimensional model for a transverse dune field which 
reproduces the observation that the dunes in a field attain approximately 
the same heights and velocities. The fact that our results do not depend on 
the initial configuration of the field and on the details of coalescence effects 
between dunes is of particular interest, since the genesis of dune fields is still 
not understood by geographers and geologists. Of course, serious limitations 
emerge from considering dunes as points in a line. Further work is required to 
compare the results presented here with real dune field observations. 
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